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ABSTRACT 

This  report  presents  an  analysis  of  radar  range  tracking  in 
the  presence  of  noise.  Three  approaches  to  the  problem  are  pre- 
sented. The  first  (Section  Tl)  uses  the  mean  square  error  crite- 
rion to  obtain  an  optimum  system.  The  second  (Section  III)  takes 
into  account  the  finite  range  aperture  of  the  system  and  an  optimum 
system  is  obtained  using,  a criterion  of  performance  which  is  more 
significant  for  range  tracking  than  the  mean  square  error  criterion. 
In  Section  IV,  a solution  is  presented  for  the  probability  of  losing 
the  target  in  the  course  of  a track.  This  last  analysis  presents 
some  mathematical  difficulties  and  seems  possible  only  in  a track- 
ing system  with  position  memory. 


RADAR  RADGE  TRACKING  AND  NOISE 


I . Automatic  Range  Tracking . 

The  basic  technique  used  for  automatic  range  tracking,  involves 
the  comparison  of  the  true  target  ranee  with  the  tracking  system 
output,  which  in  general  deviates  slirhtl.;/  from  the  true  target 
range.  The  difference  of  the  two  signals  controls  the  system  output. 
The  true  range  R is  usually  measured  by  the  time  of  return  of  a ra- 
dar echo  and  the  error  is  measured  by  the  difference  in  the  time  de- 
lay of  the  received  echo  and  a locally  generated  waveform, * 

A simple  block  diagram  for  a tracking  system  is  shown  in  Figure 
1,  The  transfer  functions  of  the  blocks  are  K and  A(s),  as  indicated 
where  ’s’  represents  the  Laplace  Transform  parameter. 
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Figure  1 

The  error  detector  has  an  output  proportional  to  the  error 
(R-R1)  with  a constant  of  proportionality  dependent  on  signal 
strength.  This  can  be  seen  by  examining  the  operation  of  a typical 
unit  which  operates  in  the  manner  described  below. 

Figure  2 shows  the  transmitted  and  received  pulses  as  well  as 
two  locally  generated  gates  which  are  delayed  in  time  from  the  trans- 
mitted pulse  by  an  amount  Rr.  For  convenience  a range  scale  rather 

*A  detailed  discussion  of  radar  tracking  systems  my  bn  found  in 
Reference  1, 
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ABSTRACT 


This  report  presents  an  analysis  of  radar  range  tracking  in 
the  presence  of  noise.  Three  approaches  to  the  problem  are  pre- 
sented. The  first  (Section  II)  uses  the  mean  square  error  crite- 
rion to  obtain  an  optimum  system.  The  second  (Section  III)  takes 
into  account  the  finite  range  aperture  of  the  system  and  an  optimum 
system  is  obtained  using  a criterion  of  performance  which  is  more 
significant  for  range  tracking  than  the  mean  square  error  criterion. 
In  Section  IV,  a solution  is  presented  for  the  probability  of  losing 
the  target  in  the  course  of  a track.  This  last  analysis  presents 
some  mathematical  difficulties  and  seems  possible  only  in  a track- 
ing system  with  position  memory. 
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RADAR  RANGE  TRACKING  AND  NOISE 
I.  Automatic  Range  Tracking. 

The  basic  technique  used  for  automatic  range  tracking  Involves 
the  comparison  of  the  true  target  range  with  the  tracking  system 
output,  which  in  general  deviates  slirhtly  from  the  true  target 
range.  The  difference  of  the  two  signals  controls  the  system  output. 
The  true  range  R is  usually  measured  by  the  time  of  return  of  a ra- 
dar echo  and  the  error  is  measured  by  the  difference  in  the  time  de- 
lay of  the  received  echo  and  a locally  generated  wavef orm.*:- 

A simple  block  diagram  for  a tracking  system  is  shown  in  Figure 
1.  The  transfer  functions  of  the  blocks  are  K and  A(s)#  as  indicated 
where  ’s*  represents  the  Laplace  Transform  parameter. 
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Figure  1 

The  error  detector  has  an  output  proportional  to  the  error 
(R-R*)  with  a constant  of  proportionality  dependent  on  signal 
strength.  This  can  be  seen  by  examining  the  operation  of  a typical 
unit  which  operates  In  the  manner  described  below. 

Figure  2 shows  the  transmitted  and  received  pulses  as  well  as 
two  locally  generated  gates  which  are  delayed  in  time  from  the  trans- 
mitted pulse  by  an  amount  R*.  For  convenience  a range  scale  rather 


*A  detailed  discussion  of  radar  tracking  system*  may  be  found  in 
Reference  1, 
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than  a time  scale  is  used  and  units  of  time  are  converted  to  their 

equivalents  in  range,  ; 
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Figure  2 

The  action  of  the  error  detector  is  to  integrate  the  received 
signal  in  the  interval  of  the  early  gate  and  then  in  the  interval  of 
the  late  gate.  The  difference  of  those  integrals  constitutes  the 
error  signal.  It  is  clear  that  the  discriminator  output  depends  on 
both  the  signal  and  error  amplitudes.  Figure  3 shows  a series  of 
error  detector  characteristics  for  different  signal  strengths.  Ex- 
perimental curves  closely  resemble  those  theoretical  characteristics 
but  do  not  exhibit  the  discontinuities  which  are  due  to  the  infinite 
slopo  of  the  leading  and  lagging  edges  of  the  pulse. 
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Error  Detector  Characteristic  (TR  = T) 
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The  value  of  K in  Figure  1 la  the  slope  of  the  oharaoterlstio 
in  the  normal  oporating  region  and  is  given  by 
(1-1)  K = 2 S fr 

where  S ia  the  signal  amplitude  whon  gated  and  fr  is  the  pulse  repe- 
tition frequency. 

The  function  unit  transfer  characteristic  A(s)  is  determined  by 
the  type  of  response  desired  and  is  almost  always  either  one  or  the 
other  of  the  following* 

a)  A(  s)  = — (an  integrator) 

Cs 

b)  a(s)  *s  — (a  single  and  double  integrator) 

Cs  CDs*- 

where  C and  D are  constants  of  the  system. 

Upon  a complete  signal  fade  and  in  the  absence  of  noise,  the 
output  of  system  (a)  will  remain  constant,  exhibiting  position  memory, 
while  the  output  of  system  (b)  will  vary  linearly  with  time.  This 
linear  variation  corresponds  to  the  velocity  of  the  target  (for  a 
constant  velocity  target)  and  in  this  manner  system  (b)  exhibits 
velocity  memory,*- 

The  overall  transfer  function  for  the  traoklng  system  is 


(1-2) 

This  loads  to 
(1-3) 


F(s)  = 


F(  3 ) = 


KA(s) 

1 + KA(s) 


ts  + 1 

for  the  position  memory  system,  and  to 


d-4) 


r,/  \ 

F(s)  = — -A 2— 

X a ♦ S + w'rx 
n 


•k 


' t 

i«4 


I 

1 

t 


♦From  this  point  on  system  (a)  will  bo  roferrod  to  as  the  Position 
Memory  or  PM  System  and  system  (b)  as  the  Velocity  Memory  or  IfM 
System. 
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for  the  velocity  memory  system,  whore  t --  — in  tho  time  constant  of 


■"  -Jh 


is  tho  undamped  natural 


tho  position  memory  system  and 
frequency  in  radians  of  tho  velocity  memory  system# 

Tho  steady  stato  position  error  for  a target  of  constant  velo- 
city ’v*  is 


v-r  for  PM  System 

o for  VM  System 

Tho  steady  state  position  error  for  a target  of  constant  accel- 
eration *af  is 


(1-5)  eov  = Llm  jf{l-F(a>\ 

S — >0  i 


(1-6) 


oa 


Lim 

s-*o 


fz((i  " F(s)}ss 


for  PM  System 
for  VM  System 


In  the  pres  once  of  noise  the.  video  echo  appears  as  in  Figure  If 

so  that  the  orror  detector  output  now  contains  in  addition  to  the 


Figure  if 

error  signal  a random  component  due  to  noise#  The  random  component 
is  independent  from  pulse  to  pulse  since  the  correlation  time  of  the 
video  noise  is  of  the  order  of  the  reciprocal  of  the  receiver  band- 
width which  is  much  loss  than  the  pulse  repetition  period.  If  v/o 
assume  that  the  moan  of  tho  gated  output  is  due  entirely  to  the 
presence  of  the  pulsed  signal  and  the  fluctuation  is  duo  entirely  to 
the  poise  so  that  signal  and  noiso  affect  tho  system  independently, 
we  may  write  for  the  output  of  tho  error  detootor 
(1-7)  1 = 1,  + 1„ 

whore  la  la  K(  R-R ' ) and  : ..  la  tho  random  component  duo  to  nolaoc  Wo 
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alao  modify  our  blook  diagram  of  Figure  1 to  that  of  Figure 

WHtev 


D 

Error 

Detector 

L 

Function 

Unit 
A(s ) 

— >■ 

K 

i=K(R-RH 

I 

Figure  5 

ijj  consists  of  a series  of  pulses  of  duration  much  smaller  than  the 
pulse  repetition  period  and  which  have  a mean  value  zero.  The  power 
spectrum  of  a pulse  series  of  this  type  has  been  found  by  Uhlenbeck# 
to  be 

(1-8)  G(f)  =2  fr  j B ( f ) | 2 (a2) 


where  fr  » pulse  repetition  frequency 

B(f)  = amplitude  spectrum  of  a single  pulse  of  unit  size 

and  a2  = mean  square  value  of  the  pulse  amplitude. 

The  spectrum  has  a width  of  the  order  of  the  reciprocal  of  the 

pul3e  duration  and  the  tracking  system  has  a pass  band  of  the  order 

of  one  cycle  so  that  for  our  purpose  we  shall  be  primarily  interested 

in  the  spectrum  in  the  region  near  f = 0.  Now  G(o)  may  be  interpreted 
2 2 

to  be  2 fr  og  where  og  is  the  variance  of  the  integral  of  the  noise 

p 

pulse,  so  that  if  iN  is  a current,  0g"  is  the  variance  of  the  random 

p 

charge  for  each  gating  operation.  The  value  of  <Jg  for  an  early- 
late  gate  system  and  square  law  envelope  detector  can  be  shown  to  ba 
approximately^:* 


*See  Section  of  Reference  2, 
aaSaa  References  3 and  If. 
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(1-9) 


T 


2 


6 


2 


6 


N 


where  d ^ is  the  dispersion  of  the  video  signal  at  the  cate.  Equa- 
tion (1-9)  is  valid  only  for  T_  rouchly  equal  to  the  reciprocal  of 

O 

the  receiver  bandwidth.  This  is  usually  the  case  in  practical  sys- 
tems , 


II , Mean  Square  Error, 

We  now  turn  to  the  calculation  of  the  fluctuation  of  the  output 
R’.  If  we  assume  that  the  system  operates  only  on  the  linear  portion 
of  the  error  detector  characteristic  then  this  calculation  may  be 
performed  by  first  finding  the  overall  transfer  characteristic  F^(s) 
of  ipj-  to  the  output  R and  from  this  obtain  the  spectral  density  of 
the  output.  The  integral  of  the  soectral  density  will  give  us  the 
output  fluctuation.  The  calculations  will  be  made  for  a VM  system 
and  the  results  for  a PM  system  can  be  obtained  from  these  by  letting 
D ■ > oo  or  <*>n ajo, 

Fjj(s)  can  be  computed  with  the  aid  of  Figure  5,  to  be 

P 


(2-1)  Fn(s)  = 

The  spectra 

is  0(f)  I Fk( J2  n 


A(s ) _ F( s ) _ _1_  3 + un  T 

1 + K A(s)  K K ts  2 + s + un2T 

? density  of  R»  (excluding  the  & -function  at  f = 
2 

f)|  , so  that  the  dispersion  of  R>  is  given  by 


0) 


.-on 


(2-2) 


O2  = 


G(0)|FH(  J2n  f)|  df  = 


2u  o 


6 


bO  .On  O O 

(4  + u»  2 tc)df 

i ^ n 


K2 


o 


kn  2?2  + (u>2n^  - 4it2f2i:  )2 


The  integral  represents  the  equivalent  noise  bandwidth  b of  the 
system. 


An  evaluation  of  the  Integral  gives 


6 


(2-3) 


b - 


1 + u) 

k' 


2 _2 


n 


2fr,  a 


The  term  r JjQ  may  be  computed  from  (1-1)  and  (1-9)  in  terms 

K2 


of  signal  and  noise  amplitudes  in  the  video  giving 


(2-10 


,2=  V 

21. 


N 


It  should  be  noted,  however,  that  d is  not  directly  proportional 
to  the  reciprocal  of  S since,  for  a given  system,  b is  also  a func- 
tion of  S. 

We  are  now  prepared  to  design  a system  to  have  minimum  error 
for  a fixed  signal  to  noise  ratio.  The  system  parameter  values  which 
give  low  noise  fluctuations  also  result  in  high  steady  state  errors 
[see  equations  (1-5),  (1-6)  and  (2-i|)]  , therefore  the  total  mean 
square  error  should  bo  considered.  This  Is  given  by 


(2-5) 


e = e, 


where  eQ  is  the  steady  state  error. 

Consider  first  a PM  system  to  be  designed  for  minimum  error  for 

a target  with  velocity  v and  a video  noise  to  signal  power  ratio  of 
0 2 

N , From  (1-5),  (2-3),  (2-i+)  and  (2-5),  we  get 

2 o 

P 2 2 ^ p1  ® m 

- Vc  1 + ■-*. u 


8fr 


Minimizing  this  with  respect  to  t we  get  for  minimum  e the  condition 
that 


(2-6) 


x = 


T 2 2 il/3 

V>  0 K 


Li6frs2v2  J 


For  this  value  of  x the  fluctuation  due  to  noise  is  twice  the  mean 
square  steady  state  error* 

The  VM  system  has  two  parameters  u>n  and  t which  may  be  adjusted 
for  minimum  error.  For  a target  with  acceleration  *a»  we  get  from 
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(1-6),  (2-3),  (20ij.)  and  (2-5)  that 

,2 

r-  -fi-  - 

S 


T 2 1 + W 2 T2 

2 _ _a  __2_  n. 

TF  2f- 


. r 

and  from  this  expression  we  obtain  the  conditions  for  minimum  £ as 
( 2-7 ) un  t:  = 1 


co  = 
n 


f 16  a2  fr  S2] 

Tr,a  a / J 


i/5 


If  is  fixed  from  other  considerations,  then  for  minimum  e with 

respect  to  variation  in  x wo  still  pet  u x = 1. 

This  represents  an  underdamped  system  [u>nT  = l/2  corresponds  to 
critical  damping  (see  equation  1-ip ) ] . Since  those  results  would 
generally  be  applied  to  systems  with  low  signal  levels  it  should  be 
noted  that  unx  varies  as  the  reciprocal  of  the  square  root  of 
signal  level  and  for  increased  signal  amplitudes  the  damping  in- 
creases . 


Ill,  Another  Crlterlor  for  System  Optimization, 

The  analysis  and  system  optimization  of  Section  II  in  no  way 
takes  into  account  the  fact  that  the  system  is  not  linear  throughout 
the  entire  range  of  error  values  (see  Fig.  3).  The  range  aperture 
for  linear  operation  extends  about  the  value  R*  = R and  is  determined 
by  the  gate  and  echo  widths,  A simple  technique  which  we  may  use  to 
take  cognizance  of  the  finiteness  of  the  range  aperture  is  to  quan- 
tize the  error  by  associating  one  value  with  it  if  it  falls  within 
the  range  aperture  and  another  If  it  is  outside  this  aperture.  This 
contrasts  to  the  mean  square  error  which  weights  the  error  values  in 
proportion  to  the  square  of  the  error  magnitude. 

Here,  as  in  the  previous  section,  we  shall  assume  complete  linear 
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operation  and  our  first  task  is  to  find  the  probability  of  the 
error  falling  within  the  range  aperture.  Since  the  noise  input 
to  the  tracking  system  is  wide  band  compared  to  the  band  of  the 
tracking  system,  the  output  R and  the  error  e will  be  Gaussian 
distributed,  e has  a mean  value  eQ  given  by  (1-5)  and  (1-6)  and 

p 

a dispersion  O'  given  by  (2-I4),  so  that  the  probability  M that  the 
error  be  within  the  ranee  aperture  is  given  by 


(3-D 


SZh  d 


+L/2 


-L/2 


(e  -e  ) 
' c' 


where  L is  the  magnitude  of  the  range  aperture  centered  at  e=  0. 

The  value  of  L may  be  made  equal  to  T (see  Fig.  3)  but 
in  an  application  for  which  It  is  important  to  restrict  the  error 
to  a specific  interval,  L may  be  made  equal  to  this  interval, 
liquation  (3-1)  ma'r  .:o  lewritten  as 

L/2-e0  _ r2 

0-2)  m = -=  r ° o dr 

-L/2-e0 


where  r is  the  normalised  variable  — — — . 

0 

For  low  signal  levels  where  fq  is  rather  close  to  L/2  and 
the  value  of  M may  be  maximized  by  maximi?ing  the  value  of 

L/2  - c 


This  is  equivalent  to  finding  the  values  of  wn  and  x for 


which 
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<3-3 ) 2 £t(eo°2)  + (L/Z  ‘ 3e<=>  3r  = 0 

and 

<3-4 ) 2 -j—  (e^2)  + <L/2  - 3 e0 ) = 0 

The  results  for  the  P.M.  system  are  obtained  by  substitut- 
ing (1-5)  (2-3)  and  (2-1+)  in  (3-3)  and  solving  for  optimin  x giving 


(3-5) 


x = 


L 


Tv~ 


The  results  for  the  V .1* • system  are  obtained  by  substitut- 
ing (1-6)  (2-3)  and  (2-1+)  in  (3-3)  and  (3-4)  and  solving  for  optimum 
u>n  and  X giving 


(3-6) 


“n't=  1 


(3-7) 


U) 


n 


Condition  (3-7)  is  valid  regardless  of  any  constraints  on 

the  value  of  u t but  the  optimum  value  of  w x is  given  by  (3-6). 
n n 


IV.  Probability  of  Target  Loss. 

The  analyses  of  Sections  II  and  III  consider  steady  state  or 
stationary  solutions  for  the  output  of  the  tracking  system.  In 
reality  a steady  state  is  never  reached  if  ore  takes  into  account 
the  true  nature  of  the  error  detector  characteristic.  For  any 
system  with  noise  the  error  will  at  times  reach  a sufficiently  large 
value  to  give  zoro  output  from  the  error  detector.  ’Then  this  occurs 

**The  result  given  by  this  equation  was  found  by  Tuchmore  et  al  (see 
reference  5)  for  a critically  dnnped  system  using  a different 
criterion  for  optimization. 
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the  existence  of  a target  echo  doos  not  act  as  r restoring  "force" 
and  the  system  output  is  governed  by  noire  only.  The  noise  may 
act  to  reduce  the  error  and  to  bring  the  system  under  control  of 
the  target  echo  again  but  It  may,  also,  act  to  increase  the  error. 
Ultimately  the  output  will  leave  the  aperture  not  to  return. 

An  important  question  that  may  be  asked  is,  "'..'hat  is  the 
probability  that  the  target  is  not  lost  in  time  t?"  By  "lost" 
we  mean  that  the  error  does  not  exceed  a certain  limit  for  the 
whole  period.  The  limit  nay  be  the  limit  of  linear  operations  or 
may  be  a value  determined  from  other  considerations  for  a particular 
application. 

This  problem  is  a rjrmidable  one  and  a solution  seems  possible 
only  in  the  case  of  r P.M.  system  with  complete  linear  operation 
within  the  limits  imposed  on  the  error. 

The  method  of  solution  i3  as  follows: 

1.  Set  up  the  Fokker-Planck”  equation  for  the  probability  density 
of  the  output  R1,  This  partial  differential  equation  is  derived 
from  the  system  transfer  function. 

2.  Solve  the  equation  with  the  boundary  conditions  which  specify 

the  output  at  the  start  of  the  track  and  the  fact  that  the  system 

is  "absorbed"  if  the  error  exceeds  the  imposed  limits. 

3.  Integrate  the  solution  for  the  probability  density  in  the 

restricted  interval  to  obtain  the  probability  of  the  system  not 
having  left  this  interval. 

The  differential  equation  which  describes  the  P.M.  system  may 
be  obtained  from  Figure  5 as 

complete  discussion  of  tho  method  of  Fokker-Planck  and  of 
Brownian  ?!otion  problems  may  be  found  In  references  6 and  7# 
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TdR'  + R«  = R + 
dt  K 


The  solution  for  R»  may  be  considered  in  two  parts--the  first  as  the 

response  to  the  forcing  function  R and  the  second  as  the  response  to 

the  noise  1^.  Let  the  first  bo  R»  and  the  second  x#  Thus  for  x we 

have  the  differential  equation 

(4-1)  dx  + x = ijj 

dt  % K 1 

This  equation  is  analogous  to  the  Lnngevin  equation  of  Brownian 
Motion, 

The  derivation  of  the  Fokkor-I-lanck  equation  from  ()_|-1)  will  not 

L 

be  given  here  but  may  be  found  in  the  articles  of  Chandrasekhar0  and 
’.Vang  and  rthlcnbeck^  . The  Fokker-Planek  equation  for  the  probability 
density  is 


(4-2) 


where  is  the  dispersion  of  R*  (and,  also  x)  found  in  Section  II  for 

a completely  linear  error  detector  characteristic.  The  term  cs 

corresponds  to  one  quarter  of  the  spectral  density  of  ill  which  is  the 

Kt 

appropriate  coefficiont  for  the  second  derivative  term  In  (4-2) 
according  to  the  method  of  Fokkor-Planck, 

The  first  thing  wo  can  do  is  to  check  our  results  for  the  case 
handled  In  Section  II,  This  case  corresponds  to  imposing  the  follow- 
ing boundary  conditions  on  Ui -2): 

(4-3)  p (±  «*,  t)  = o 

p (x,o  ) * 5(x  - x ) 

o 

The  solution  with  those  bounds ry  conditions  is  found  in  Appendix 
A and  is  given  by 

(MU  p(x,t>  = [2*02(l-O  -2‘A)]-i/2  oxpr.  (x-x0a  1 

\ 2s  ) 


This  corresponds  to  a Gaussian  distribution  for  x v'.th  a tino 
varying  moan  given  by  XGe  "^A  and  n time  varying  dispersion  given  by 
<j2(l  -e  The  steady  state  vain,  a ire  and  0^  respectively 

which  agree  with  results  in  Section  II, 

It  is  Interesting  to  note  what  haupens  to  p(x,t)  for  the  limiting 
cf.se  of  zero  noise . A Gaussian  probability  density  with  its  dispersion 
approaching  zero  is  in  the  limiting  case  a delta  function  so  that 

Urn  ? p(x,t)  = &(  x - xQe  *fcA) 

d^o 

This  tells  us  in  effect  that  if  the  noise  should  suddenly  cease 
at  t = o and  at  that  time  x = xQ  , thereafter  x ’"ill  have  a value 
xoe  “tA  with  certainty.  This  is  in  exact  agreement  with  our  know- 
ledge of  the  normal  transient  response  for  the  system, 

'Vo  now  turn  to  the  problem  of  the  solution  for  the  probability 
of  not  losing  the  target  in  the  sense  discussed  at  the  beginning  of 
this  section,  '• 'e  shall  do  this  for  a zero  velocity  target  so  that 
eo  = o and  x - e , For  this  we  must  find  a method  of  eliminating 
from  our  probability  calculations  the  possibility  that  8 = R -R1  bo 
outside  the  range  aperture  at  any  time  in  the  tricking  interval,  Ve 
shall  consider  a target  lost  if  at  my  time  in  the  observation  time 
interval,  the  noise  is  sufficient  to  cause  a deviation  of  the  error 
e such  that  |e|  >,  L/2, 

To  calculate  the  probability  of  losing  a target,  we  may  employ 
an  artifice  in  the  problem.  Let  us  consider  the  case  where  two 
"absorbing"  vails  are  placed  at  either  end  of  the  range  aperture  (soo 
Figure  6),  These  walls  have  the  characteristic  that  if  R*  should 
have  a value  equal  to  R = L/2  tho  system  would  * i absorbed  and  not 
released.  Therefore,  with  the  absorbing,  walls,  If,  after  some  tine, 
e is  found  within  the  rango  aperture,  then  it  necessarily  remained 
within  the  rango  aperture  the  entire  interval j otherwise,  it  would 
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! ) 

i 


l6 


t 


have  been  absorbed  and  "out  of  circulation".  The  existence  of  the 
walls  has  no  effect  on  tho  fluctuations  of  R»  as  lone  as  e remains 
within  the  aperture.  The  probability  of  not  losing  tho  target  in 
our  problem  is,  therefore,  tho  same  as  that  of  R not  being  absorbed 
by  one  of  the  walls.  Tho  effect  of  the  absorbing  walls  is  to  intro- 
duce tho  boundary  conditions  that  for  xQ  = o, the  probability  density 
p (x,t)  = o for  x « - L/2.  Mathematically,  the  problem  remains  to 
solve  with  the  boundary  conditions  that 

P {•  h , t)  = 0 

p ( X , o ) = b ( X ) 


The  solution  for  low  signal  levels  is  given  in  Appendix  B and 
the  result  is  that  tho  probability  of  not  losing  the  target  is 
a-S)  P=£  ±-  e-"2e*!i|Qnm 


n=l ,3,5  n* 


n 


whero  ^ and  0 are  dimensionless  parameters  givon  by 


p p 

® ~ ft  0 t and  it  « L 

L2  x F3 


and  Qn  ) is  defined  by 

°n  <♦)  « 


nn 

r'~r  . + 5 

e n*-  it c cos  x dx 


Tho  function  P ( ^ ,9)  is  tabulated  as  a function  of  9 for 
several  different  values  of  in  Appendix  D and  is  grep  hed  in 
Figure  7, 

For  higher  signal  levels  a different  method  of  solving  (1-2)  is 
necessary.  This  is  done  in  Aopendix  C.  The  results  are  tabulated  in 
Appendix  D and  graphed  in  Figure  8. 


) 

l 
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Appendix  A Solution  of  the  Fokker-Plonck  Equation  in  an  Infinite 
Interval. 

The  boundary  conditions  are  that  p (-°o  , fc)  = 0 and  that  at 
t = 0,  x = xQ  with  certainty.  This  is  equivalent  to  p(x,0)  » 6(x 
Tho  solution  may  be  found  as  follows} 

Change  variables  of  (l;-2)  to 

z = x e^/T 
t»  = t 

Since 

JSL  = -IE  + JL.  et//x  "q| 

8t  fit*  i 


and 


= e 

O 


2tA  92p 

az 


the  original  equation,  after  dropping  the  prime  from  t becomes 


ft  P = 1 p + _C_ 

at  t -T 


2 2t/T  JLR 


ft  2 


with  the  initial  condition  that  at  t = 0 

p =s  ft(z-x0) 

Now,  if  wo  let 

p = efc/T  f(z,t) 

we  obtain  for  f(z,t)  the  equation, 


e2t/x  2%  ; fU.oj  = 6U-*0) 


ftt 


We  may  now  separate  variables  by  letting  f(z,t)  * S(z)  T(t) 
with  tho  result  that 


II  0-2tA  3 sn 


T a' 


- A , a constant 
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Prom  this  wo  pot 
T = (constant)  exp 

and 


o 2 u2tA  dt  i = ( constant  )oxp(-X2d^(  o2t^« 

_ J t T 


0 1 


S = constant  o 


j\z 


Since  the  constants  in  the  solution  for  S and  T depend  on  X and 
all  values  of  X are  permitted,  wo  may  write 

f(z,t)  = 1 A (X ) oxp  J jXz  -X2  (e2tA  - 1) 

V — C5o  2 

From  tho  condition  that  f(z,0)  = ft(z-x^)  wo  have 


dX 


,+°o 

A(X  ) e 


JXz 


dX  =6  (z-x0) 


or 


A (X)  = 


;“JXxc 
2 it” 


This  value  for  A(X  ) may  bo  substituted  in  tho  expression  for  f(z,t) 
and  after  integration  the  result  is 

i2  i 

\ 2d2(02tA  -i)  / 

Making  the  necessary  substitutions,  wo  obtain  for  p as  a function 
of  x and  t the  expression 


f(z,t)  = [2*o2  (G2tA:  -I)j  ~l/z  cp,  - ( |~xg  j- 


p(x,t)  = [aw^l-o-2*/*  ) ] -V2  oxp  ! - -~-0-  .a'2t/^  )2 1 

j 2o2  (1-o-2‘/t,  / 
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- Solution  of  Fokkor-Planck  Equation  In  a Finite  Interval 


and  Low  Signal  Levels, 

The  boundary  conditions  on  (Ii-2)  are  given  by 

P (tL/  ,t)  = 0 

p (x,o)  = 6 (x) 

The  coefficient  1 depends  only  on  signal  strength  and  not  at  all  on 

X 

noise,  while  c£  depends  only  on  the  noise  level  and  not  on  signal. 

X 

For  zero  signal  \ vanishes  and  ( ip—  6 ) with  its  boundary  conditions 

x 

reduces  to  a standard  problem  in  linear  heat  conduction,  (e.g.  see 
Churchill8,  p.  102) 

•Vo  may  separate  variables  by  writing  p (x,t)  in  the  form 
(B-l ) p(x,t)  = v(x)‘  u(t) 

so  that  for  v(x)  and  u(t)  v/c  obtain 

-Z  — = iw  *■  + — = -X  , a constant, 

o ^ u <J^  v v 


- Xe  ^ t 


From  this  we  got 

\ ft  2 4- 

(B-2)  u **  constant  o ” ° - 

x 

and  for  v(x) 

(B-3) 

v”  + ill  + (i_  + X)  v = 0 
0 c 0d 

The  boundary  conditions  on(ij-2)  imply  that  v (t  =0 

For  low  signal  levels  1 is  small  and  we  may  solve  (B-2)  to 

T? 

a first  order  in  1 • But  to  do  this,  it  is  necessary  to  put  (B-3) 

in  a form  from  which  it  is  possible  to  eliminate  second  order  terms* 
This  may  bo  dono  by  using  a new  dependent  variable  f(x)  such  that 

(b-4)  . . . =£  . . 


v(x)  » o'  IJJ*  f (x) 


ao  that 


v"(x) 


o"  5e2 
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Now  (B-3)  In  terms  of  f(x)  bocomes 


(B-5) 


with  the  boundary  condition 

f (1  L/2)  = 0 

v2f 

Let  us,  for  the  time  being,  ignore  the  1 , term  and  we  shall 

4 ■ 

examine  hater  under  what  conditions  this  is  warranted.  Equation 
(B-5)  now  remains  as 


(B-6) 


f"  + (-L-  + X)  f = 0 ; f (t  h = 0 

2?  2 


The  solution  is 


<B-7)  i = « ooa  x | i + x.  + 3 sin  x /~i — r~jT 

Vso5  -Jio5 

where  A and  3 are  constants.  The  condition  f (i  ^)  = 0 

2 

requires  that  B = 0 and  that 


f - A cos  x / l + \ + 3 sin  x /“T  “ 


(B-8) 


I 1 + X = -D.7L-  ; n = 1,  3,  5,  

[2o2  L 


Therefore,  X.  may  only  assume  the  values 

(B-9)  X = ; n * 1,  3,  5 

L2  2c- 

These  represent  the  eigenvalues  of  X for  equation  (4-2)  with  its 
associated  boundary  conditions. 

The  comolete  solution  may  be  found  by  combining  (B«4),  (B-7), 
(B-8)  with  (4-2)  to  obtain 


(B-10) 


p(x,t) 


^ 7 f « x1-  _ n2i i 2 o2  t + t ) 

Z_  Ane*p  — FT-  t Rf 

n=l,3,5-..  1 ) 


It  nay  be  noted  that  the  lowest  eigenvalue  X^  is  equal  to 

2 

■■  *■  • 1 and  the  highest  value  of  x in  (B-5)  occurs  at  the  boundary 

l>  2c* 

2 2 2 

x « L/2  , Therefore,  in  (B-5),  — Xr_  l>  and  1 , + \ v *_ 

WT'iSfr  Sc? 
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so  that  neglecting  th©  former  term  does  not  introduce  much  error 


if  Lc 


« * 


The  constant  An  can  be  determined  from  the  initial  condition 


that  p(x,o)  = 6(x)  or 


(B-U) 


OQ  ^ 

An  e"  coa  ILL*  = &(x) 


n = 1,3,5 


_ xJ 


Multiplying  both  aides  by  e [^2  coa  m 2L  - (m  = odd) 

L 

integrating  from  - L/2  to  L/2  , and  recognizing  that 

-L/2  L/2  n = m i 

COS  m 7'  x COS  n^JTCJK  dX  as  V 

il/2  L L ° n * m ) 


cos  m 7'  x cos  n 7i  x dx  = 
L L 


we  have 


A = £ 


To  find  the  probability,  P,  of  not  losing  the  target  we  must 
integrate  the  probability  density  p (x,t)  in  the  range  aperture, 
(B-32)  *L/2  _ L/2  p 


p(x,t )dx  = 


-L/2 


, ^ L/2  2 

n=l,3,5..  -L/2 


We  define  the  dimensionlcaa  quantities 

0 = x2  £ and 


L*  * 


(B-13) 


\|>  = JL 
2 c 


which,  when  substituted  in  (B-12)  gives 
(B-lW  ....  «=“-  l 


n-1,3,5 


-wThe  function  P(o#«®  ) has  been  tabulated  by  Olson  and  Schultz 
in  Reference  9*  ** 


i 
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where  the  function  Qn(iJ> ) is  defined  by 

( B-i5)  JL2L  4,2  x2 

Qn(^)  = f 2 e cos  x <jx 

J o 

The  function  Qn(\J> ) is  briefly  discussed  and  tabulated  in  Appendix  D. 


Appendix  C - Solution  of  Fokkor-Planck  Equation  in  a Finite  Interval 
with  Higher  Signal  Levels, 

Hero,  as  in  Appondix  B,  the  boundary  conditions  on  (ij.-2)  are 
given  by 

p(t  L/2 , t ) = 0 
p(x,o)  - 6 (x) 


and  as  before  vve 

(B-2)  and  (B-5). 
(C-l) 


start  by  separating  variables  and  obtaining  equations 
In  equation  (B-5),  substitute 


y = 


x 

o 


In  terms  of  the  now  variable  y,  wo  obtain 
( c-2) 

1 * X ’ )f  » o 


d2f 


+ 


where 


X 1 = X 


with  the  boundary  conditions  that  at  y = ± |l.  s + 4,  f f = 0 

Equation  (C-2)  Is  the  Weber  equation  and  can  be  found  discussed 

in  the  literature,  (e,g,  see  Ince  10  , Weber11  , Whittaker  and 

12 

Watson  ,)  The  solution  to  this  equation  for  integral  values  of  \» 
is  the  parabolic  cylinder  function  vhich  may  be  expressed  in  closed 
form  as  -1/2 

(c-3)  = tx'  : 2 rr*i 


-^A 


o 


Vy) 


2$ 


where  It  . »(y)  la  the  Hormitc  polynomial  of  degree 
\ 

For  non-integral  values  of  \ wo  may  solve  equation  (C-2)  but 
not  in  closed  form  by  first  defining  a now  dependont  variable  by 

r2. 


(C-4) 


f = e 


-y 


^4-  w(y) 


thereby  transforming,  equation  (C-2)  to 


(c-5) 


d2w 


dy 


2 


We  may  now  assume  a series  solution  for  w, 

(C-6)  W = 1 + h^y  + l^y2  + h^y3  + . . . . 

Substituting  this  in  (C-5),  we  find  that  the  coefficient  of  yn 

is  (n  + 2)(n  + 1)  h ^ 0 - nh_  + \'h  . Setting  all  the  coefficients 

n + c.  n n 

equal  to  zero,  we  obtain  a recurrence  formula 


(C-7) 


_ X 


*n  + 2 


+ n 


(n+1) (n+2 ) 

This  leads  to  two  solutions  for  w--cne,  a power  series  containing  all 
terms  even  powered  in  y,  and  one  with  all  odd  powered  terms  in  y, 
giving  an  even  function  end  odd  function  solution  respectively . 

Since  our  result  is  necessarily  symmetrical  about  the  zero  axis,  wo 
are  left  with  the  solution 

(c"8>  W , X . X'  2 + \'U'  - 2)  y4  . V(X'-2)(\'-lt)  6 

v U:  6; 

It  now  remains  to  find  the  values  of  X*  for  which  w vanishes  at  the 
boundaries  y = ± ^ • A method  applicable  hero  was  used  by  Chand- 
rasekhar1^- in  a problem  in  Stellar  Dynamics*  Ho  attributes  this 
method  to  Sonmorfold1^* 


oThe  functions  of  the  parabolic  oylinder  may  bo  found  plotted  in 
Jahnko  and  Erade1^  p,33. 
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If  wo  oxami.no  (C-2)  with  reference  to  our  discussion  of  the 
parent  oquation  (B-0),  wo  may  note  that  the  solution  for  (C-2)  for 
small  ij)  can  be  accomplished  by  neglecting  the  term  y . This 

leads  to  a sorios  of  eigenvalues  for  the  lowest  of  vhich  is 


(C-9) 


\ = 
l 


n 


TT2  ' 2 


This  lowest  eigenvaluo  decreases  with  increasing  \|>  • For  t|)  , 

wo  cannot  use  the  value  in  (C-9)  but  may  obtain  the  eigenvalues  with 
same  reasoning  used  often  in  Quantum  Mechanics  to  find  eigenvalues, 
(o.g.,  see  the  discussion  of  the  harmonic  oscillator  in  Pauling  and 

Wilson^,  p.  71),  The  reasoning  is  as  follows:  The  solution  f(y) 

- iv2 

in  equation  (C-2)  is  e Ir  w(y)  where  w(y)  is  given  by  the  power 

series  ^C-8).  This  result  will  not  converge  for  y — ^.ooin  spite  of 

the  e ^ factor,  unless  the  series  ( C— 8 ) terminates  after  a finite 

number  of  terms.  This  may  be  demonstrated  by  comparing  the  series 

y 2/2 

for  w(y)  to  the  series  for  e , 


(C-10)  y2/ 2 


‘*r*sS;* 


The  ratio  of  the  (n  + 2 )—  term  coefficient  to  the 
for  large  n is  ^ , 

h +2 

The  recursion  formula  (C-7)  also  gives  as  the  ratio  -2 — for 

1 hn 
large  n the  value  — , Therefore,  coefficients  for  the  higher  degree 

n 

terms  of  y in  (C-8)  and  (C-10)  differ  only  by  a constant  ratio,  so 
that  w(y)  is  of  the  order  of  e^  ^for  large  y*  Since  f (y)*e~^ ^ 
f(y)  will  diverge  for  y — 

For  convergence,  (C-8)  must  terminate  after  a finite  number  of 

terras,  or  X * must  take  on  the  values  0,  2,  if These  represent 

the  eigenvalues  of  X ' for  ^ lowest  eigenvalue  being  sero« 


th 


term  coefficient 
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It  would  acom  therefore,  that  the  lowest  eigenvalue  takes  on  tho 
values  given  in  (C-9)  for  low  values  of  ^ and  approaches  0 as  t|» — 
as  shown  in  the  sketch  on  page  28.  The  second  lowest  eigenvalue  by 
the  same  token  will  behave  similarly  and  approach  2 as  an  asymptotic 
value , 

V/ith  this  In  mind,  let  us  rearrange  the  terms  of  (0-8)  so  that 
our  solution  is  in  terms  of  a power  series  in  X*  , We  have  then 
(C-ll)  w(y ) = i -X’q^Cy)  + \'  2 <p2(y)  - X*2  <p2(y)  - ^^(y)  . . 


The  functions 
(C-12) 


<P1(y),  <P2(y)  ..  are  power  series  in  y: 

<P1(y)  = a2y2  + a^y^  + a^y6  + ... 

+ b6y^  + 


etc . 


Since  w = 0 at  y =\|i  we  may  substitute  this  condition  in  (C-ll) 
and  solve  for  w(y)  by  first  considering  only  the  first  two  terms  in 
the  series  and  neglecting  the  higher  order  terms  in  , obtaining 
for  the  lowest  eigenvalue 
(0-13) 

1 *,.<♦> 


Greater  accuracy  may  bo  obtained  by  including  another  term  in  the 
series  (C-ll)  and  obtaining  ^ 

(C-ll*)  X 91(^}  13(^}  ' ^92U)1 

1 2 92(  ) 

For  values  of  ^greater  than  2,  ( C — ) and  even  (013)  offer 
very  good  accuracy.  The  results  of  these  calculations  are  tabulated 
In  Appendix  D,  For  small,  the  method  breaks  down  because  it  would 
be  necessary  to  work  with  many  terms  of  equation  (C-ll).  This 
presents  diffieulties  in  tho  computation  of  the  terms  and  then  in  tho 
solution  of  a high  degree  polynomial  in  It  is  interesting  to  note 
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that  for  4' » 1,  tho  lowest  eigonvaluc  X*  , by  the  method  outlined 

above  gives 

X ^ - 2 • 03^4-  using  equation  (C-l4)»  and 

« 1*999  by  Including  thoX'  term  in  (C-ll)  and  solving 

. » 

the  resulting  cubic  equation  in  , 

i 2 , 

while  by  the  method  of  Appendix  B,  we  obtain  X = — - — = 1.967 

1 4 2 

This  excellent  agreement  indicates  an  overlapping  in  range  of 
applicability  of  the  two  methods  of  solution,  thus  making  unnecessary 
any  higher  order  approximation,  such  as  Perturbation  or  W.K.B.  methods* 
to  extend  the  range  of  the  method  of  solution  used  for  low  ^ in 
Appendix  B. 

To  each  eigenvalue  x’n  corresponds  an  eigenfunction  Wn(y)  , 
which  is  a solution  to  (C-5)  and  an  eigen  function  f (y)  = e ^^%n(y), 
which  is  a solution  to  (C-2).  Tho  homogeneous  differential  equation 
(C-2)  is  in  standard  Sturm-Liouville  form  with  homogeneous  boundary 
conditions  so  that  the  functions  fn(y)  form  an  orthogonal  set. 

°r,  p*  p*  _y2/2 

(C-15)  l fn(y)  fn(y)  dya  I e Wn(y)  wm(y)  dy  » 0 n / m 

-t|>  U 

The  solution  to  the  original  partial  differential  equation  may 
now  be  formed  by  combining  (B-2)  , (B 4*.)  * (C-l)  and  (C-4-)  into  (B— 1 ) • 
The  solution  is 

(°-l6)  p<X,t>  = C A 9 'K’n  t/Z  O*  k2  Wn  (*) 

n n n d' 

The  constant  ^ may  be  determined  from  tho  initial  condition 

p(x,o)  * 6 (x)  or 

( C-17 ) r — > jjL 

^Ane  5*  »„<£>“»<*> 

®The  proof  of  this  statement*  as  well  as  adiacussion  on  Sturm-Liouvllle 
problems  may  be  found  in  Churchill®,  Ince*”,  f'argenau  and  Murphy^®* etc* 
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In  terras  of  y = ^ , this  is 

<c-l8)  CA 


V”"-  = i »(y) 


Since  fn(y)  = e 


-y% 


Wn(y)  form  an  orthogonal  set,  if  we  multiply 


In  other  words, 

* 


both  sides  of  (018)  by  Wm(y)  and  integrate  over  the  interval  y = 

to  y = +\[i  , all  the  terms  in  the  summation  will  vanish  except  that 
for  n = m 

n -L 

■+ 

leads  to 

(c-19)  V 1 


p +tji  n +4* 

n • Wn(y)  Wm(y)dy  * I &(y)  Wm(y)  dy 

H.J. 


2 o 


w^(y)dy 


V/e  can  find  the  probability  of  not  losing  the  target  from  (C-16) 
and  (C-19)  by 

TL/2 

P * 1 p(x, t ) dx 

-L4 


(C-20) 


Or  in  terms  of  the  variable  y 
( C-2X ) _ f*+  v . » 


\ 


n 


A_  o"  Xn  * e"'  ''Z  w„(y)  Odj 


n 


-tl» 


Substituting  the  value  of  An  leads  to 

(C-22)  -.-XA  Vt  it  e‘y2,6  Wn(y)  dy 

0 


£ 

n 


f*  .-»%■§<  y> 

i!  o 


uy 


Since  the  first  eigenvalue  Is  near  zero  aid  the  second  in  the 
order  of  two,  etc.,  we  see  that  all  terms  except  the  first  in  equation 
(C-22)  will  diminish  rapidly  to  zero  with  increasing  t/x  • 
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Ample  accuracy  ia  obtained  if  we  use  only  the  first  term  of  the 
aeries,  or  ^ 

<o-23)  . x!  v»  1 n°'y2/2  Wl(y)dy  _ ,-x: 


P * e 


r e“y2/£  W*  (y )dy 


* Cx(  tl»)  e"*l  ^ 


A tabulation  of  the  coefficient,  C^(i|>  ) , as  a function  of  \|>  may 


be  found  in  Appendix  D, 

It  may  seem  surprising  that  the  coefficient  has  a value  very 
slightly  greater  than  one,  which  would  seem  to  indicate  that  for 
0 — ; P > 1 , a situation  which  is  impossible.  The  difficulty  lies 
in  the  fact  that  only  the  first  term  in  equation  (C -22)  was  used. 

The  second  eigen  function  would  give  a negative  contribution.  At  any 
rate,  this  introduces  an  error  only  for  V*  small  a*  discussed  above. 


Appendix  D 

1,  Evaluation  ofQn(\}i) 

The  function  0 D 

mt 

GnW>)  = e”  oos  x dx  n = 1,3,5 

nay  be  found  in  terms  of  Error  Integrals  with  complex  arguments. 
Complete  tables  for  this  are  not  available  and  therefore  the  function 
had  to  be  evaluated.  Since  our  Interest  is  primarily  In  low  values 
of  4 it  was  found  advisable  to  expand  the  exponential  function  as  a 
power  series  in  \|>  and  evaluate  terra  by  term,  each  term  being  of  the 

* fnV2 

*0rm  Jo  X 003  x d*  (n  °dd).  This  integral  is  tabulated  for  n 
oven  in  reference  18,  but  for  n odd  a brief  tabulation  sufficient  for 
our  purposes  is  shown  below. 
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* 

n = 1 

ro 

II 

C 

1 

n = 5 

0.1 

.9995 

-.9977 

.9976 

0.2 

.9981 

-.9909 

.99014- 

o.5 

.9983 

-.9M|.i 

.94H 

1.0 

.9550 

-.7876 

.7819 

2.  Tabulation  of  P (\|> , 0 ) from  (4-5) 


A. 

9 

PW»  . 9) 

e 

= 1.0 

PW>.9) 

0 

1.000 

0 

1.000 

0.2 

.975 

0.2 

.981 

o.L 

.8  pc 

0.4 

.874 

0.6 

.710 

0.6 

.752 

o.e 

.588 

.487 

0.8 

.642 

1.0 

1.0 

.548 

1.5 

.303 

1.5 

.368 

2.0 

.166 

2.0 

.247 

Tabulation 

of  Eigenvalue 

Xj  as  a function 

of  ( see 

Appendix 

ij)  First  order  approx.  Second  order  apnrox. 


2 

.222143 

.22 685 

6 

.088013 

.09277 

3 

*023491 

.02395 

10 

.015105 

.01531 

4 

.000994 

.000996 

Evaluation  and 

The  coefficient 

Tabulation  of  (See 

C^(i|0  is  determined  as 

WjCyMy- 

J 0 

Appendix  C) 

^ -y^/2  2 ' 

6 2 V (y)dy 

o 
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whero  w-^(y) 
term  in  the 
We  may  find 


is  a power  series  in  y and  X.  ’ , Therefore,  a 

r+ 

evaluation  involvos  an  integral  of  the  form  3 

J o 

tabulated  in  the  literature*  the  function 


!!n(x)  * 


12  tt(n-i):: 


y11  0-^2  dy 


from  which  a typical  term  in  the  evaluation  of  c (i|i)can 


c1(t|>  ) 


be 


2 1.0835 

6 1*0538 

3 1.0208 


10  1.0160 

I4.  1.0011 


typical 
ne'y2/2  dy. 


computed. 


g~ - " ■ . . 

The  function  ^(x)  for  values  of  n below  11  is  tabulated  in 
reference  19,  and  for  n » 11,12  in  reference  20. 
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